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Abstract 

We present a simple model for fermion mass matrices and quark mixing in the context of 
supersymmetric grand unified theories and show its agreement with experiment. Our model 
realizes the GUT mass relations = 3me, rris = mi, = irtj. in a new way and is easily 

consistent with values of rrit suggested by MSSM fits to LEP data. 
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Understanding the values of fermion masses and quark mixing remains one of the most 
important unsolved problems in particle physics. The magnitude of the challenge is shown 
by the fact that at present it is not definitely clear which approach one should take to 
the related matter of electroweak symmetry breaking (EWSB). Intensive efforts continue in 
at least two quite different directions: (i) supersymmetry, and, in particular, the minimal 
supersymmetric standard model (MSSM), in which EWSB is perturbative (ii) dynamical 
(and nonperturbative) EWSB, e.g. slowly running technicolor 0, as well as others such as 
compositeness. Here we shall adopt the first approach, specifically the MSSM, which may 
provide the simplest way to understand the precision electroweak measurements from LEP |Q 
while stabilizing the Higgs sector in the standard model. A further motivation for the MSSM 
is that it achieves unification of gauge couplings at a single scale rriG, a prerequisite for a 
supersymmetric grand unified theory (SGUT). 

In GUT's such as SU(5) or SO(IO), by well-known restrictions on (renormalizable) Higgs 
couplings, the elements of the lepton Yukawa matrix can be made equal to 1 or —3 times 
corresponding elements of down quark Yukawa matrix. In this approach, one can get a 
simple relation between the resultant (running) lepton and down quark masses, evaluated at 
the GUT scale. One interesting proposal is that at the GUT scale 



It is important to determine which (experimentally acceptable) forms for fermion mass ma- 
trices yield this relation. Since 1979, essentially only one such form (up to trivial equivalence 
transformations) has been found. With -C^n = Vui^u,Lyu'^u,R+Vde{'4'd,Lyd'^d,R+'4'e,Lyei^e,R) + 
h.c, this form has GUT scale Yukawa matrices 
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Interesting studies of SGUT Yukawa matrices, and this ansatz in particular, have been carried 
out recently using renormalization group equations (RGE's) in the MSSM |^, (where 
Vu = Vde = 2-i/2ycos/5, with v = = 246 GeV). A general feature 

of this ansatz is that it tends to require a rather large top quark mass, rrit. Further recent 
works have considered more complicated mass matrices P] and higher-dimension operators |Q 
which do not, in general, yield (|ip (see also Ref. |T^). 

In this paper we exhibit a new and quite different model for fermion mass matrices 
which yields the GJ mass relation (|l]). We show that this model agrees with experiment, in 
particular with values of nit around 135 GeV, in the range suggested by MSSM fits to the 
ni. Our model is defined at the SGUT scale by 
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A natural setting for this ansatz would be an SO(IO) SGUT, where one can obtain symmetric 
(complex) Yukawa matrices by using 10 and 126 Higgs representations to couple to the 
fermion 16 x 16 bilinear. If one had a fundamental theory of fermion masses, one would 
presumably be able to derive the forms of the Yf and also the values of the parameters from 
first principles. In the absence of such a theory, we believe that models such as ours can 
give valuable hints about the underlying physics []12[. Our mass matrix model has seven real 
parameters Au, B^, Cu, Ad, Bd, Cd, and (p (plus tan/5) to describe the nine fermion masses in 
the u, d, e sectors and the four angles parametrizing the quark mixing matrix V (plus tan /3); 
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hence it makes six predictions. (Because of the uncertainties regarding neutrino masses and 
mixing, we do not consider these here.) 

We first show that this ansatz yields the GUT scale GJ mass relation (|I]). The Yf, f = 
u, d, e are diagonalized by unitary transformations Ufi UfYfUj = diag{\f^i, where 
A/j = ±mfj/vf (with the quantities defined at /i = mc, and Vd = Ve = Vde)- Without loss of 
generality, one may pick > 0. For one can take either of the choices (A^^i, \d,2^ A^^s) = 
(+,-,+), (-,+,+); correspondingly, (Ae,i, Ae,2, A^^s) = (-,+,+), (+,-,+). To leading 
order in small fermion mass ratios, we find that Cd = Ce ^ mh = rrir. Next, Bd = 
A(i,2 = T^s/vde and Ae,2 = — 3Ad^2, which together imply = Finally, Ad = 



Ad,iAd,2 = ^rndrris/vde and (^6)12 = {Ye)2i = ^-K,iK,2 = ^/mjri^/vde so that the 
property (Ve)i2 = Ad implies ^Jm^fn^ = ^Jnidml and hence rud = Snie. (These relations 
hold for either set of sign choices.) This proves our assertion. We will use the choice 
(— , +, +) for Yd since it gives a slightly better fit to the data on quark mixing. For Yu one 
may also take (A„,i, A„,2, A„,3) = (+,-,+) or (-,+,+), and one finds C„ = A„,3 = rrit/vu, 



Bu = A„,2 = ^fnc/vu, and Au = y —Xu,i^u,2 = y/fnum-c/vu- The different sign choices for the 
Xuj will imply picking values in the fit which differ by vr; for definiteness, we use (+,—,+) 
for the Xuj- 

To test the model, we evolve the Yukawa matrices down from the SGUT scale to the 
electroweak scale, ttLj^^ and diagonalize them there. We shall take rrij^^ = rrit, but our 
results are not very sensitive to this choice. We use the illustrative value mt{mt) = 135 GeV 
(i.e., pole mass Mt = mt{mt)[l + Aas^Mt) / (Sn) + 0{a3{Mt)'^] ^ 142 GeV), consistent with 
MSSM fits to LEP data [0, To avoid complications involving further model-dependence, 
we follow a common simplification |^ of approximating the full SUSY mass spectrum by 
a single mass scale, nisusY with msusY = ^t- Given the theoretical uncertainties 
associated with the SUSY spectrum and threshold corrections, 1-loop RGE's are sufficiently 
accurate in the range from to msusY- The inputs Ciemi^z)^^ = 127.9±0.2 and sin^ 9]y = 
0.2325 ± 0.008 [|14|, which give ai(mz)"^ = 58.90 ± 0.11 and a2("^z)~^ = 29.74 ± 0.11, 
together with the assumptions that the lightest Higgs has rrih — mz and msusv — ^t = 135 
GeV lead to 1-loop unification in the MSSM at mc = 1.3 x 10^^ with a^^ = 24.8. We 
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take a^{mz) = 0.125, consistent with current measurements 0, |T^] and with MSSM gauge 



couphng unification. (Using the 3- loop QCD RGB's JT6|, this gives = 302 MeV and 



as{mt) = 0.118.) 

An interesting property of the model which we have established by analytic and numerical 
solutions of the RGB's is that although the values of the parameters change, the forms of 
Yu and Y^. are preserved to high accuracy under evolution from to rrij^^. For Y^, the 
relation (Yd) 21 = (^d)i2 the equality (Yd) 22 = (Xd)23 are maintained quite accurately, 
while 1(1^)32! decreases (by about 35 %) relative to | (1^)23!. Keeping dominant terms (and 
with f(t) = (1671^)^^(1/ /dt, t = ln(/i/mG)), the RGB's for the Yukawa matrices thus reduce 
to RGB's for the nonzero elements, which are 



for Xu = Au, Bu] 



for Xd = Ad, (Fd)22, (Ya 



d)23, 



for Zd = (Yd) 32, (Yd) 33; and 



X^ = Xy(-J2cy,jg]) (8) 
j 

Cu = Cy(-J2cy,g^ + QCl) (9) 
j 

Xd = Xd(-J2^d,j9]) (10) 
j 

Zd = Zd(-j:^d,,9] + Cl) (11) 
j 

A:, = Xe(-Ece,,^7|) (12) 



for Xe = Ae, Be, Ce, whcrc for the MSSM in the relevant energy range /i > msusY, 

(6l,62,&3) = (f , 1, -3), (Cu,l,Cu,2,Cu,3) = (^,^,y)^ (Cd,l> Crf,2, £^,3) = (^, 3, y ), (Ce,l, Ce,2, Ce.s) 
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3, 0) |jl^|. To this order, the running of is negligibly small. 



Manipulating these RGB's, we obtain the solutions 

mb(m^^) Cd(t^^-) fCu(tEw)V'^ 



^Amjsw) Ce(t„^) \ Cu(tG) 
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andt^,^ = \n{mt/mG). We calculate (t^,^) = 1.92. Further, we define 77/ = mf{mf) /mfinij^ 
for leptons and heavy quarks, and rjuds = ""^/(l GeV)/mf{mj^^) for u, d, s and take the value 
mb{mb) = 4.19 GeV, consistent with Ref. 0. Using 1-loop QED and 3-loop QCD RGB's, we 
obtain rjr = 1.02 and rjb = 1.55, whence mb^rrij^^) /mr^rrij^^) = 1.54. We successfully fit this 
ratio by choosing an appropriate value of Cuitc)- This is a nontrivial test of the ansatz, since 
it is not, a priori, guaranteed that a perturbatively acceptable value of Cuitc) would enable 
one to fit the observed mb/rrir mass ratio. We find Cu{tj^y^.)/Cu{tG) = 0.269 from (|T3|); substi- 
tuting this into the solution to the RGE for Cu [0, we get Cuitc) = 4.48. (This is consistent 
with the validity of perturbation theory since loop corrections are small: C^/(167r^) = 0.13.) 
We then get Cu{tJ^y^,-) = 1.21 and, from the condition mt{mt) = Cu(t^^y)f sin iP)/V2, deter- 
mine tan (3 = 0.84. 

From the RGB's, we obtain the solutions 

Ad{t) Bait) 
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3 

,=\v«,(t) 

We calculate pAsit = t,^^) = 2.38. Now 

) Bdit 

EW 



pABit) (15) 



pAB{t) = m:^) " (16) 



0.792 (17) 



^l^i^Ew) 3i?e(t 

EW 

Using m^(m^^y) = 102.8 MeV, the model predicts nisirrLj^^) = 81.4 MeV, i.e., with our 
3-loop rjuds = 2.94, 771^(1 GeV) = 239 MeV. This is consistent with the (upper range of the) 
determination ms(l GeV) = 175 ± 55 MeV |2^. Further we find 



(>d)23(t) \Cuitc)J fiAa,(t) 



3 ".j 
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We use the results Adit)vde = \Jmdit)msit), Aeit)vde = y'me(t)m^(t) and eq. (P!3|) to find 
^dim^^-)/meim^^) = 7.13. Using me(m^^y) = 0.4876 MeV, this yields mdirrij^^) = 3.48 
MeV, and hence mrf(l GeV) = 10.2 MeV, in agreement with the value mrf(l GeV) = 8.9 ±2.6 
MeV \M. 
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Given that Yu retains its form to high accuracy under the evolution from tug to rrij^^, 
the two remaining parameters in Yu are determined via the relations Bu{tew)vu = —mc{t) 



and Au(tew)vu = \/ ^u(t)'mc(t) . From mdrnc) = 1.27±0.05 GeV we calculate rjc = 2.32, 
whence m^[m^^) = 0.547±0.022 GeV. With m„(l GeV) = 5.1±1.5 MeV we get m„(m^^) = 
1.7 ± 0.5 MeV. Diagonalizing at m^^, we have vJJu,lYuUIj^ = M„, vlUd,LYdY}ulj^ = M|, 
whence V = Uu,lU\i for the Cabibbo-Kobayashi-Maskawa (CKM) matrix V PT|. To leading 
order \Vus\ — \\l'mu/^c — e^'^^-^md/msl (with quantities evaluated at = m^^^ [EH)- 
determine (p = 75.5° by fitting \Vus\- Quoting as usual the rephasing- invariant \Vjk\ values, 
we get 

0.9753 0.221 0.002 
{\Vij\} = I 0.221 0.9748 0.030 | (19) 
0.006 0.030 0.9995 



in agreement with current data [T^ (with \Vcb\ in the lower range of acceptable values). 
For the reparametrization-invariant CP violation parameter J |j2^, we get J = 0.95 x 10~^, 



again in agreement with data. Further details will be given elsewhere [^. 
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